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Abstract 



The behavior of a two dimensional, steady turbulent mixing layer was investigated. Be- 
sides the usual velocity components, all the contributing components of the Reynolds 
Stresses are also calculated and presented. The results indicated that in the two dimen- 
^ 'sional steady turbulent mixing layers the theory of self-preservation is not valid for all the 
flow domain, and that the flow is turbulent regime is divided into two areas. Through 
O .calculation of Reynolds Stress components at the point of initiation of the flow, it was 
^ shown that the turbulent flow in the area containing the point of singularity of the flow 
•^h 'is not, as previously believed, self-preserved. 
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g 1 Introduction 

Turbulence is the most frequently occurring mode of fluid motion. It is characterized by 
high fluid particle mixing, and high energy dissipation [1-3]. The exact time-dependent 
nature of events leading to the onset and development of turbulence is not fully known; 
however, by the aid of time averaging techniques we can obtain meaningful and practically 
useful statistical definitions [3]. Nearly all the problems that are studied in the field of 
Turbulent Fluid Dynamics are numerically extensive and expensive [3-5], and there is 
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no single concept to model different classes of turbulent flow regimes [3-5]. In general a 
turbulent flow is characterized by random variation of flow variables, high levels of fluid 
particle mixing, and high energy dissipation [2, 3]. 

Turbulent flow regimes may be divided into two main categories: Free Turbulence, 
where no physical boundary limits the development of flow: Bounded Turbulence, where 
flow is restricted by physical boundaries [2]. Experimentally, Free turbulent flows are 
easy to generate, and that is primarily why they have been so extensively studied. From 
a theoretical stand point, their formulation and simulation is also much simpler, since 
for normal incompressible flows, no pressure gradient term has relevance [3]. In studying 
free turbulent flow fields it is desirable for practical reasons to confine oneself, at first, 
to rather simple but important elementary flow patterns [1, 2]. These elementary flow 
patterns may be classified into two main groups: (1) free jets, and (2) wake flows behind 
obstacles. The flow pattern under investigation here is a plane, two dimensional mixing 
layer [3]. In this flow pattern it is possible to distinguish one main flow direction with 
velocity substantially greater than in any other direction [1, 2]. 

Through out the past few decades, turbulent plane mixing layer has been the subject 
of numerous experimental and theoretical research, and it is widely used to investigate 
various aspects of phenomena related to free turbulence, including process and nature of 
mixing [3], structure of turbulent components and their self-preservation [3-6], stability 
of the flow [4], transfer of momentum and energy [5], intensity, and rates of temporal and 
spatial expansions of free turbulent flows [6]. 

The spatial and temporal development of the turbulence is closely correlated to the 
development of Reynolds Stress Components. The existing literature indicates that for 
free turbulent flows, the entire flow domain was considered to be self-preserved [7]. How- 
ever, in this investigation it is shown that the self-preservation theory is not applicable to 
all the flow field, and that at areas very close to the point of initiation of turbulent flow 
regime the said theory is not applicable, i.e. the flow is not self-preserved. 

Following the introduction the fluid dynamics equations are indicated and the formu- 
lation is completed. In section three, the simulation results and discussions are presented. 
And in the final section discussion and concluding remarks are presented. 



2 Turbulent plane mixing layer's governing equations 

For a 2-dimensional incompressible turbulent plane mixing layer, the governing equations 
are Reynolds equations and the continuity equation. They may be written in the following 
form: 
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where bar represents time-averaged quantities, and u±, u 2 are the velocity components in 
the x and y directions respectively; also u[, and u' 2 are the fluctuating components of the 
Ui and u 2 velocity components respectively. 

To close the above set of equations the Reynolds Stress Equation Model ("RSM") 
is used [4]. The exact equation for the transport of Reynolds Stresses, R^, takes the 
following form: 

OR 

-^ + diviRijV) = Pij + Dij - dj + lUj + tin , (4) 

where P^, D^, e^, 11^, and represent rate of production of Rij, transport of Rij 
by diffusion, rate of dissipation of R^, transport of R^ due to turbulent pressure-strain 
interactions, and transport of R^ due to rotation respectively These terms are represented 
by the following equations: 

p..--(R. ^L + R- ^) (5) 

where V bar indicates the velocity vector, and having 

„ d ,u t dRij . . . 

Dij = t— (— o-^) , (6) 

OX m Ufa ox m 

also, 

v t = {k 2 /e)C, , (7) 

in which = 0.09, and = 1.00. 

Assuming isotropy of the small dissipative eddies, and assuming that it affects the 
normal Reynolds stresses (i = j) only, and in equal measures we have [4]: 

en = (2/3)6^ , (8) 

where the Kronecker delta, 5ij is given by 5ij = 1 if % = j and Sij — if % ^ j. 

The transport of Reynolds stresses due to pressure-strain interactions is calculated by 
the following expression: 

n« = -Ci(|)(fltf - (2/3)foS«) - C 2 (P« - (2/3)P%) , (9) 
in which C\ = 1.80 and C 2 = 0.60. The rotational term is given by: 

2ix>fc (Rjrrfcikm Rim^jkm) ■ (10) 

Here Uk is the rotation vector and is the alternating symbol; = 1 if i, j and k are 
different and in cyclic order, = — 1 if i, j and k are different and in anti-cyclic order 
and e^k = if any two indices are the same. 

Turbulent kinetic energy k is needed in the above formula and can be found by adding 
the three normal stresses together: 

k = 0.5(R U + R 22 + R 33 ) (11) 
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and the equation for the transport of scalar dissipation rate e is: 

^ + div(peV) = div(^grade) + 2(e/k)C U fi t E ij ■ E tJ - p(e 2 /k)C 2£ , (12) 
at <y £ 

where C± £ = 1.44 and C 2e = 1.92, and Eij is the mean rate of deformation of a fluid 
element and E^- • Eij is their scalar product [4] . Note that for the free stream velocities 
Rij = and e — 0. 

Here the values of the above-defined relationships are presented. For R 12 we have: 

p (r dU2 + R dUl Mn dU2 + R dUl \ nr> 

fi2 — — h-Ki2 7i r-«42T; r n 2 2T. ) 
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also, 
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n 12 = (-1.8e/fc)i2i 2 - 0.6P 12 (18) 
fii2 = (19) 
£12 = 0. (20) 



For Ru components we have: 



P U = .2(^ + ^*1), ( 21 ) 

OX 1 OX2 

Dn = W^'W^ + ^ 2 { To^ h (22) 

£11 = (2/3)e, (23) 

n n = (-1.8e/ife)(i2ii-(2/3)A;)-0.6Pii (24) 

fin = 0. (25) 



And for R 22 components we have: 



P22 = -2(^22^ + ^2^), (26) 
OTi OX 2 
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d u t dR 22 d v t dR 22 

° 22 - ^ ( io^r )+ ^ ( LoW' (27) 

e 22 = (2/3)e, (28) 

n 22 = (-1.8e/fc)(i22 2 -(2/3)A;)-0.6P22, (29) 

fi 22 = 0. (30) 

For £ we have the following expressions: 

de drj de d . v t de . d , u t de . 



where 



= ^, (32) 



B 22 = (33) 
dx 2 

El2 = 0.5(^1 + ^), (34) 



and the kinetic energy of turbulence is 

k = 0.5(^11 + #22) . (35) 

Therefore we have obtained six equations which must be solved simultaneously to 
provide the two unknown velocity components, three Reynolds stress components and 
the scalar turbulent dissipation rate. 

With the aid of standard Computational Fluid Dynamic (CFD) schemes and formu- 
lations, namely the Central Finite Difference Discretization, and using a forward time- 
marching scheme, the equations are numerically simulated and solved. The time marching 
is carried out until the steady state conditions were established. The x, y and t steps are 
determined by trial and error and are fixed and equal to 0.0001/40 m, 0.00001/40 m, and 
10~ 7 seconds respectively. 



3 Simulation results and discussions 

The respective computer codes for the CFD formulation of the turbulent plane mixing 
layer are run until the steady state solution is obtained. Using the approximation, for 
free stream velocities of U max = 10.0(m/s), U min = 5.0(m/s), the velocity component 
distributions and the non-dimensional Reynolds stress components at various flow sec- 
tions are obtained, and for sections located at about x = O.OOOOl(m), near the initiation 
point of the flow, where inherently a critical point of discontinuity exists, and about 
x = 0.000036(m), far down-stream from the point of initiation of the flow, are plotted 
and presented in the following figures. In the following figures U re f = U max — U min . 
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Figure 1 indicates the experimental results for the spatial expansion of the turbulent 
mixing layer obtained by D. Oster and I. Wygnanski [7]. In the same figure, the simulation 
results obtained from the utilized CFD code are indicated. 

In Figure 2, the Reynolds Stress distribution obtained experimentally [7], at x — 
300mm, x = 500 mm through 1700 mm, and by the utilized CFD formulation, at the 
area close to the point of initiation of the turbulent flow, is shown. 

Figures 3 and 4 indicate the variation of u velocity component against the non- dimen- 
sional vertical distance at two intersections close to the point of initiation of the mixing 
layer, as calculated by the utilized CFD scheme. 

Figures 5 and 6 show the distribution of non-dimensional u' 2 Reynolds Stress Compo- 
nents at cross sections close to the point of initiation of turbulent flow, calculated by the 
used CFD scheme. 

From these figures, the reduction in the peak value of turbulent intensity in the flow 
direction indicates that the in the region under consideration the flow is not self- preserv- 
ing. 

Figures 7 and 8 present the variation of the non-dimensional v 12 Reynolds Stress Com- 
ponents at cross sections close to the point of initiation of turbulent flow, determined by 
the used CFD code. 

Figures 7 and 8 also indicate that the turbulent flow at the desired area is not self- 
preserving. In addition, comparing the order of magnitudes of u' 2 and v 12 Reynolds Stress 
components in figures 8 through 10, although the relative contribution of u' 2 component 
to the generation of turbulence is more, its effect on further development of turbulence 
regime downstream of the point of initiation of the plane mixing layer is reduced. 

Figures 9 and 10 show the distribution of the time averaged, non-dimensional auto- 
correlated components of the turbulent velocity fluctuations, very close to the point of 
emergence of the mixing layer, obtained by the CFD solution. 

From figures 9 and 10, comparing the maximum value of the turbulence intensity with 
the experimental values measured by Oster and Wygnanski [7], Yule, and Spencer, both 
presented at the same reference and indicated in Table 1, the acceptable performance of 
the CFD solution and our results is established. 

Table 1. The experimental and theoretical maximum value of intensity 
for U max /U min ratio of 0.5. 



Source 


u'v'/U 2 ef 


Spencer (1970) 


0.011 


Yule (1971) 


0.013 


Oster and Wygnanski (1982) 


0.013 


Present result 


0.013 



In addition, Figures 5 through 10 indicate that turbulence Reynolds Stresses, which 
are themselves responsible for turbulence generation, have their peak value at the center 
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of symmetry of the flow. This fact together with the order of magnitude of Reynolds 
Stress variation complies with the known behavior of such turbulent flow regimes [3-7]. 

Based upon their experimental set up, erected to study free and forced turbulent plane 
mixing layers, D. Oster and I. Wygnanski had deducted that for free stream velocity ratios 
of less than 0.6 the flow remains self-preserving [7]. However, their area of interest started 
200 mm down stream of the point of initiation of the turbulent flow. 

The results obtained by the utilized formulation indicated that for the area of interest 
in such flow patterns the auto-correlated Reynolds Stress components have closer spatial 
behavior, but the cross-correlation of turbulent velocity components has a different con- 
tribution to the generation of turbulence, as it shows a larger peak at 0.00001 m, and 
decays faster at 0.000036 m from the point of initiation of turbulent flow regime. 

Referring to the our results, D. Oster and I. Wygnanski's deduction is invalid for the 
area of the flow close to the point of initiation of turbulence regime. In addition, fig- 
ures 5 through 10 indicate that at points further away from, but close to, the point of 
initiation of turbulence flow, the auto-correlated components of Reynolds Stresses play 
the dominant role in generating and maintaining the turbulence behavior. As a further 
expansion to such deduction it can be noted that at initial stages of such flow regimes, 
turbulence starts and expands more due to contribution of cross-correlated component 
of velocity fluctuations, and as the influenced region of turbulent flow increases, most of 
the generated momentum and energy of turbulent components needed to generate and 
maintain the turbulent flow pattern is supported through the influence of the auto- corre- 
lated fluctuating velocity components. In other words, there are two regions of Reynolds 
Stress contribution to the generation and expansion of turbulence. In the first region, 
the turbulence is generated through a higher correlation of vertical fluctuating velocity 
components; such close correlation decays yielding a second mechanism of turbulence gen- 
eration, forming a region at which turbulence is mainly maintained by higher temporal 
auto-correlation of fluctuating velocity components. 

4 Conclusion 

The particular area under consideration includes a point of discontinuity, critical point, 
where all the flow variables are constantly zero [3-7]. As the starting point of the turbulent 
flow regime, the consequent development of the turbulent mixing layer is dependent to 
events that take place in this very small region. Comparing the results with similar 
experimental research on free turbulent plane mixing layers, we deduct that the obtained 
results are acceptable [5, 7]. However, the complex nature of events in this region and 
their ultimate effect on the development of turbulent flow deserve much more theoretical 
and experimental investigation. 
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Figure captions 



Fig. 1. Comparison of the spreading of unforced mixing layer obtained by the utilized 
CFD and experimental set up of D. Oster and I. Wygnanski. 



Fig. 2. Distribution of time averaged \u'v'\/U^ for U max /U min of less than 0.6, provided 
by D. Oster and I. Wygnanski, and calculated by the utilized CFD scheme. 

Fig. 3. Graph of Velocity Component u, as calculated by the CFD scheme for flow cross- 
section at x = 0.00001 m, for the free stream velocity ratio of 1/2. 

Fig. 4. Graph of velocity component u, as calculated by the CFD scheme for flow cross 
section at x — 0.000036 m, for the free stream velocity ratio of 1/2. 

Fig. 5. Graph of variation of average \u'u'\/U% e f, against vertical distance, as calculated 
by CFD scheme, for flow cross section at x — 0.00001 m, for free stream velocity 
ratio of 1/2. 

Fig. 6. Graph of variation of average \u'u'\/U^ e f, against vertical distance, as calculated 
by CFD scheme, for flow cross section at x — 0.000036 m, for free stream velocity 
ratio of 1/2. 

Fig. 7. Graph of variation of average \v'v'\/U^ e f, against vertical distance, as calculated 
by CFD scheme, for flow cross section at x — 0.00001 m, for free stream velocity 
ratio of 1/2. 

Fig. 8. Graph of variation of average \v'v'\/U^ e f, against vertical distance, as calculated 
by CFD scheme, for flow cross section at x — 0.000036 m, for free stream velocity 
ratio of 1/2. 

Fig. 9. Graph of variation of average \u'v'\/U? e f, against vertical distance, as calculated 
by CFD scheme, for flow cross section at x — 0.00001 m, for free stream velocity 
ratio of 1/2. 

Fig. 10. Graph of variation of average \u'v'\/U^ e f, against vertical distance, as calculated 
by CFD scheme, for flow cross section at x — 0.000036 m, for free stream velocity 
ratio of 1/2. 



9 



This figure "Figuresland2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0301066vl 



This figure "Figures3and4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0301066vl 



This figure "Figures5and6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0301066vl 



This figure "Figures7and8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0301066vl 



This figure "Figures9andl0.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/physics/0301066vl 



